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ABSTRACT 

In previous calculations for the luminosities of 22 GHz water masers, the pumping is reduced and ultimately 
quenched with increasing depth into the gas because of trapping of the infrared («30-15() /xm), spectral line 
radiation of the water molecule. When the absorption (and reemission) of infrared radiation by dust grains is 
included, we demonstrate that the pumping is no longer quenched but remains constant with increasing optical 
depth. A temperature difference between the grains and the gas is required. Such conditions are expected to 
occur, for example, in the circumnuclear masing environments created by X-rays in active galaxies. Here, the 
calculated 22 GHz maser luminosities are increased by more than an order of magnitude. Application to the 
well-studied, circumnuclear masing disk in the galaxy NGC 4258 yields a maser luminosity near that inferred from 
observations if the observed X-ray flux is assumed to be incident onto only the inner surface of the disk. 

Subject headings: circumstellar matter — galaxies: active — masers 


1. INTRODUCTION 

Observations are accumulating to indicate that 22 GHz 
water masers frequently occur in the gas near the nuclei of 
active galaxies (e.g., Braatz, Wilson, & Henkel 1994). These 
masers have been particularly well studied in NGC 4258 
(M106) since their detection (Claussen, Heiligman, & Lo 
1984) and provide a detailed description of a thin, circum- 
nuclear disk in Keplerian motion at a radius of 0.1 pc about a 
central mass of 4 X I0 7 A/. (Watson & Wallin 1994; Miyoshi 
et al. 1995). 

In contrast to water masers in our own galaxy (e.g., Watson 
1994), the small deviations from Keplerian for the velocities of 
the individual components of the masing gas in NGC 4258 
(<1() km s l ) suggest that shocks are not the main cause for 
circumnuclear masers. Calculations have also demonstrated 
that the observed flux of X-rays can create regions of warm 
molecular gas at high densities in the circumnuclear environ- 
ments of active galaxies, as is required for water masers 
(Krolik & Lepp 1989). Neufeld, Maloney, & Conger (1994; 
hereafter NMC) have shown in detail that the necessary 
conditions for masing of the 22 GHz transition can be 
achieved and have calculated maser luminosities. The obser- 
vation that most of the maser luminosity comes from the inner 
edge of the disk in NGC 4258 is also consistent with the idea 
that X-rays directly from a compact, nuclear source create the 
masing environment. Regardless of how the environment is 
created, inelastic collisions with H 2 molecules and radiative 
cascades involving the rotational energy levels of the ground 
vibrational state (de Jong 1973) lead to the inverted popula- 
tion for the maser. The optical depths of the infrared transi- 
tions involved in the cascades are large because of 
reabsorption by other H.O molecules. With increasing optical 
depth, trapping of this radiation reduces and ultimately 
quenches the inversion. In the circumnuclear environments 
created by X-rays, this trapping significantly reduces the maser 
luminosities in the available calculations (see Fig. 2 of NMC). 
We reason in § 3 that previously calculated maser luminosities 
seem to be insufficient to create the observed luminosities, at 
least for NGC 4258, where the refined data make possible a 
close comparison. 


In this Letter we demonstrate how dust grains absorb the 
infrared radiation involved in the pumping cascade and pre- 
vent the quenching of the inversion described above. Absorp- 
tion of the infrared radiation by dust in the pumping of H : 0 
masers has previously been examined in the context of galactic 
masers (Deguchi 1981). For galactic water masers it has not 
seemed to be important because of the lower column densities 
for the regions of warm, masing gas. We focus on the 
circumnuclear disk of NGC 4258 because the detailed obser- 
vational description of these masers provides the opportunity 
for a close comparison. 

2. CALCULATIONS 

To illustrate the influence of the grains, we present a direct 
comparison with the previous calculations. We thus use meth- 
ods similar to those of NMC to calculate the emissivity <F sal 
(photons s“ ! cm 3 ) of a saturated maser as a function of depth 
into a thick slab of gas that is infinite in the directions parallel 
to the planar interface. Escape of the infrared radiation 
involved in determining the populations of the masing states is 
only allowed to occur through the interface and is treated 
using the “escape probability” approximation, in which there is 
no velocity gradient. To perform these computations, we have 
modified the code described elsewhere (Anderson & Watson 
1990) to incorporate the emission and absorption of infrared 
radiation by dust grains, to utilize static (rather than large 
velocity gradient [LVG]) escape probabilities, and to utilize 
more recent cross sections (Green, Maluendes, & McLean 
1993). By straightforward analysis, we find that the influence of 
the dust grains can be incorporated by modifying the usual rate 
equations (e.g., Anderson & Watson 1993, eq. [11]) for the 
populations in the escape probability approximation to read 

X (g,C„nj - g,C,,n,) + x + Pji) + (My - n, 

j*i )>• 

~ + /?„) + («, - >h A'#,,] = 0, (1) 

j<i 

in which n ( and n } refer to the populations (per magnetic 
substate) of states i and j of H.O in the gas, g, and g s to the 
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statistical weights, A }1 to the Einstein ,4-values, and C, y to the 
rates of excitations (or de-excitations) due to collisions with 
molecular hydrogen (see also Hummer & Rybicki 1985; 
Deguchi 1981). The summations are over states j. The terms 
that involve the escape probabilities ft (and ft) in equation 
(1) describe the emission and absorption in the absence of dust 
grains. The terms that involve ft (and /3 ! y/ ) are modifications 
caused by the presence of grains. Note that ft = ft , /3 ■ = ft, 
and S-j = S']. When the grains are much cooler than the gas, 
the source function due to grains, ft = [exp (hi\j/kT d ) - 1] 1 
goes to zero, and the terms in equation (1) that involve ft (and 
ft) disappear. 

For convenience and consistency with previous work, we 
utilize the analytic approximations of Hollenbach & McKee 
(1979; their eqs. [5.10] and [5.21]; note, however, that the 
reemission by dust was not treated there) for j3 , y and ft: 


where 


ft = {1 + t (> [2tt In (2.13 + r 2 )]" 2 } 


ft = ft + 


(1 - ft)e~ r ' 

1 + (1 - 6*)iJ 


ft> 


( 2 ) 

( 3 ) 


8 n = 


2 tty 

1 + 2a, 


In 2.72 + 


( 4 ) 


in which r i} is the optical depth to the surface of the slab at the 
center of the ij spectral line of H 2 0, r d is the optical thickness 
to the edge of the slab due to dust at the wavelength of the ij 
line, and a tj = tJt^ For i d we use an approximation based on 
Laor & Draine (1993; see their Fig. 6) as a benchmark: 
r d = t () JV(H)/ 10 23 cm" 2 , A < 50 /in i and r d = t 0 (50 gm/ 
A)tV(H)/ 10 21 cm' 2 , A > 50 pm, in which 7V(H) is the column 
density of hydrogen nuclei to the edge of the slab and t 0 = 1. 
Calculations are also performed with different r n and wave- 
length dependences to verify that modest uncertainty in the 
assumed properties of the grains leads to only modest changes 
in our results. There is, of course, uncertainty in r {/ . To obtain 
<F sar , the radiative terms in the rate equations (1) that repre- 
sent the 22 GHz maser transition are replaced by ~R(n u - nj) 
in the equation for the upper state (i = u) and by R(n lt - n,) in 
the equation for the lower state ( i = /). Equations (1) are then 
solved numerically as a function of R to obtain the limiting 
constant value <I> sat = R(n u — nj) as R — > x. In previous inves- 
tigations, <F sat has been obtained by first computing (n u - n ( ) 
from the rate equations in the unsaturated limit ( R = 0) and 
then multiplying this quantity by an effective loss rate. Poten- 
tial ambiguities arise (at least to us) in determining an effective 
loss rate. However, with the loss rate as specified in Anderson 
& Watson (1993), we find that the two methods yield similar 
<t> sai values for the regions of most interest here. 

To perform the calculations, the gas density, kinetic tem- 
perature, and H 2 0 fractional abundance within the slab of gas 
must be specified. Fortunately, these quantities have already 
been computed by NMC (see their Fig. 2) for an incident 
X-ray flux that is essentially identical with that inferred for the 
X-ray luminosity of NGC 4258 (~10 41 ergs s 1 above 1 keV; 
Makishima et al. 1994) incident onto a gas with a pressure of 
10' 1 K cm' 1 at the distance of the inner radius of the masing 
disk (0.1 pc). In Figure 1, we present the results of our 
calculations performed in the manner described above. For 
comparison, we also show the results of similar calculations 
from Figure 2 of NMC. The general agreement between our 



Fig. 1. — The log of the maser photon emissivity vs. distance into a slab of 
gas that is irradiated by an X-ray flux similar to that expected for the 
circumnuclear masing gas in the galaxy NGC 4258 at a gas (pressure/A: ) = 10* 1 
K cm \ The expected enhancement due to dust is best seen by comparing the 
solid line (NMC), which was obtained previously in a similar calculation but 
without dust, to the solid line labeled “cold dust.” Calculations are also 
presented for various temperatures T d for the dust ( T K = gas temperature). 
Our curve {dot-dashed line) best reproduces the calculations of NMC. 


calculations for with no dust and the analogous calcula- 
tions by NMC is evidence of the general consistency of the 
calculational methods (but see below). 

The influence of dust on <F sat is illustrated by our other 
calculations presented in Figure 1. It is seen most clearly by 
comparing the two solid lines which, when integrated to an 
appropriate depth to obtain the maser luminosity, differ by a 
factor of about 20 (see § 3). Calculations are also presented in 
which the dust temperature is 300 K and in which it is equal to 
the the gas kinetic temperature T K and differs from T K by 10 K. 
Modest changes (up to a factor of 3 or so) in the velocity 
dispersion Au or the dust parameter t„ will result in compa- 
rable fractional changes in <F sa( . In § 3 we reason that the grain 
temperature should be about 200 K in the masing gas for NGC 
4258, and thus the presence of dust can cause a major increase 
in the total maser luminosity. The result in Figure 1 that 
is essentially constant (except near the boundary) when dust is 
present reflects the fact that absorption is more important than 
escape for infrared spectral line radiation. Once the column 
density of H 2 0 to the boundary is large enough to cause 
sufficient scatterings that the accumulated path length of a 
photon becomes equal to the mean free path for absorption by 
dust (which occurs when r d « 0.1), further increase in the 
column density of H 2 0 (and hence the depth into the slab) is 
unimportant. After this number of scatterings, the photon is 
destroyed through absorption by dust grains. The result that 
<F sat becomes constant can also be seen in equations (1), (3), 
and (4). With increasing depth into the slab, ft dominates ft 
and becomes independent of depth, 

ft - (2rjr t} ){\n [2.72 + (t,/^ 2 

since the ratio t d /r lf is independent of depth. Equation (1) and 
its solutions must then also become independent of depth into 
the slab, as seen in Figure 1. The presence of dust thus results 
in a lower limit to the escape probability. In detail, our 
calculations for <F sat in Figure 1 are somewhat different from 
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those of NMC. This is partly compensated by our use of 
Av = 2 km s ' versus their Av = 1 km s' 1 . After the original 
version of this Letter was submitted, D. Neufeld and P. 
Maloney kindly made available the exact variation with depth 
of the physical conditions used for calculating <F sat . When these 
are used, we have verified that the effect of dust is changed 
negligibly from that indicated in Figure 1. 

3. DISCUSSION 

Most of the radiation from a compact nuclear source 
(optical, ultraviolet, soft X-rays) will be absorbed near the 
inner interface when the slab of gas has an optical depth 
greater than a few at 50 ix m, as is relevant here. This energy 
will be processed into infrared radiation that will be absorbed 
(and reemitted) to maintain the temperatures of the dust 
grains in the slab. The benchmark temperature for the dust 
grains is then the effective blackbody temperature of the 
bolometric flux F h of radiant energy incident onto the inner 
surface (= [F h /4a] l/4 ). Calculations seem to indicate that 
grains in thick slabs have essentially this temperature beyond 
an optical depth of a few tenths (see Fig. 14 of Laor & Draine 
1993). Making the assumption that the ratio of the bolometric 
luminosity to the X-ray luminosity is about 10, as is typical for 
active nuclei, then gives a dust temperature T d « 200 K within 
the masing gas at a distance of 0.13 pc from the compact 
source in NGC 4258. In contrast, when this X-ray flux is 
incident onto the interface of the slab of gas used in Figure 1, 
it creates the gas kinetic temperatures of « 400-600 K that 
are used in the calculations for Figure 1. Although the grains 
and the gas exchange energy through collisions and through 
the emission and absorption of infrared radiation, the addi- 
tional energy input to the gas due to photoionization by the 
X-rays causes the calculated temperature difference. Based on 
the above information, this difference in temperature between 
the grains and the gas is large enough that the <f> sat values 
calculated in the “cold dust' 7 limit are adequate approxima- 
tions. 

To obtain the true (i.e., integrated over angles of emission) 
maser luminosity to compare with observations, <f> sat must be 
integrated over the masing volume. For a ring of thickness h 
with an inner radius /?, the maser luminosity is 

L m = 2ir\\vRh j 9> sat dR, (5) 

where the radial extent of the masing gas is recognized to be 
much less than R. Beyond the depth at which the optical depth 
for Thomson scattering is about 1 [A(H) = 1.5 X 10 24 H nuclei 
cm -2 ], the X-ray heating is significantly decreased. We thus 
terminate the integral in equation (5) at this depth, which 
corresponds to 4 X 10 15 cm in Figure 1. For the “cold dust” 
curve in Figure 1, h v f <f> sat dR is thus 3800 L 0 pc 2 . Then, from 
the published observations of the disk in NGC 4258 [R = 0.13 
pc and (h/R) < 0.02], L m < 8 L«, in agreement with the 
observed luminosity of 7 L 0 based on the inferred beaming 
angle (Miyoshi et al. 1995, Table 1). Our value for hi/ / <f> sat dR , 
and hence for L,„ as well, exceeds that of NMC (= 180 L 
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pc -2 ), which was computed for similar conditions, by a factor 
of 20. From a similar analysis, NMC concluded that their lower 
value of hi/ / dR is compatible with the observations. The 
discrepency arises in part because the beaming angle adopted 
by NMC (11 = 27 t h/R) is smaller than that inferred from the 
observations of NGC 4258. In the interpretation of Miyoshi et 
al. (1995), the beaming angle is determined by the size of a 
nuclear continuum source, which is being amplified by the 
central masers. This source is larger than the thickness of the 
disk. Since the continuum source has not yet been detected, 
potential uncertainty remains about this interpretation for the 
beaming angle. However, H « 2irh/R is based on the premise 
that the gas that creates most of the observed maser flux (the 
nearly radial central emission) is distributed roughly uniformly 
over a radial extent of the disk that is comparable to R. In 
contrast, Miyoshi et al. (1995) seem to find that the radial 
extent of this masing gas is much smaller than R. 

Because of the narrowness of the observed spectral lines of 
the circumnuclear water masers (as small as a few kilometers 
per second), the calculations to date have only utilized simi- 
larly small velocity dispersions in calculating escape probabil- 
ities. If velocity fields with larger dispersion are present that 
are not reflected in the spectral line profiles, the calculated 
maser luminosity will be increased. Adopting a larger velocity 
dispersion in directions orthogonal to the line of sight is 
evidently one such possibility. The thinness of the observed 
disk and the apparent need for hydrostatic equilibrium does, 
however, limit such additional velocity dispersion. The fore- 
going analysis is based on the premise (as was NMC) that a 
significant X-ray flux is incident only on the inner surface of 
the disk of thickness h. If X-rays incident on the upper and 
lower surfaces of the disk due to scattering by the halo or to 
incidence at grazing angles directly from the central source are 
important, it is evident that the masing volume will be 
increased over that included in equation (5) and the maser 
luminosity will be appropriately increased. However, quanti- 
tative assessments of these possibilities have not yet been 
presented. At gas pressures greater than the 10’ 1 K cm ■ 
considered in Figure 1, the extent (in column density of 
hydrogen) of the molecular gas created by X-rays that is warm 
enough for masing is reduced. According to NMC, at (pres- 
sure/k) = I0 12 K cm 1 the gas is cool enough that masing is 
unimportant at depths into the slab for which absorption by 
dust is important in the pumping. Finally, the absorption and 
emission by dust at the wavelengths of the spectral lines of the 
gas can influence the cooling of the gas, as well as the pumping 
of the maser. However, examination of the cooling rates (Fig. 
6 of Neufeld & Kaufman 1993) suggests that this will be of 
minor importance for the conditions on which Figure 1 is 
based. 
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